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Abstract 

m . 

Cf~) , A self-consistent electronic structure calculation based on the Luttinger-Kohn model is performed 

on GaMnAs/GaAs multilayers. The Diluted Magnetic Semiconductor layers are assumed to be 
metallic and ferromagnetic. The high Mn concentration (considered as 5% in our calculation) 
makes it possible to assume the density of magnetic moments as a continuous distribution, when 
treating the magnetic interaction between holes and the localized moment on the Mn ++ sites. 
Our calculation shows the distribution of heavy holes and light holes in the structure. A strong 
spin-polarization is observed, and the charge is concentrated mostly on the GaMnAs layers, due 



to heavy and light holes with their total angular momentum aligned anti-parallel to the average 
magnetization. The charge and spin distributions are analyzed in terms of their dependence on 
the number of multilayers, the widths of the GaMnAs and GaAs layers, and the width of lateral 
GaAs layers at the borders of the structure. 

PACS numbers: 73.21.-b, 75.75.+a, 85.75.-d,72.25.Dc 
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I. INTRODUCTION 



Recent advances on the physics and technology of GaAs-based nanostructures with di- 
luted magnetic semiconductors (DMS) open a wide range of potential applications of these 
systems in integrated magneto-optoelectronic devices.— In Gai_ x Mn x As alloys substitutional 
Mn acts as an acceptor (it binds one hole), and at the same time it carries a localized mag- 
netic moment, due to its five electrons in the 3 d shell. For x near 0.05, the alloy is a metallic 
ferromagnet,— the Curie- Weiss temperature after annealing is 160K,— and the free hole con- 
centration is near 10 20-21 cm~ 3 . The possibility of application in Spintronics and Photonics 
exists because such a layer provides the injection of spin-polarized carriers in an otherwise 
non-magnetic semiconductor region of the device, eliminating the needs of a strong exter- 
nal magnetic field. Obviously, the higher the transition temperature obtained for a DMS 
thin layer, the higher the possibilities for a device with such a layer to operate near room 
temperatures. The ferromagnetic order in the metallic phase is understood, at present, as 
resulting from the indirect exchange between the Mn 2+ ions due to the local spin polariza- 
tion of the hole gas. This explanation implies the spin coherence length to be larger than the 
average distance of the localized moments. Although most of the scientific theoretical effort 
in this problem has been directed to understand the origin of the ferromagnetism, there are 
still fundamental issues to be considered in the electronic structure. For instance, the role 
played by the light holes and holes of the split-off band is a point yet to be understood. 
Such a nomenclature in GaAs is specific for bulk systems, where the tetragonal symmetry is 
preserved. In the case of heterostructures, where such symmetry is broken by the presence 
of the interfaces, a mixing occurs in the composition of the hole states. Differently from 
heavy holes, light hole states are not spin-eigenstates. In consequence, the occurrence of a 
local magnetic field is also a factor that contributes to the mixing in the composition of the 
hole states. This can be seen easily in the framework of the effective mass approximation 
and Luttinger-Kohn (LK) k-p expansion,— because off-diagonal terms appear in the Hamil- 
tonian, breaking the tetragonal symmetry. Therefore, the presence of interfaces, together 
with a local magnetic field, claims for a better treatment of the calculation of the electronic 
properties in DMS heterostructures. 

To the present, the six bands k-p method has been used to obtain the valence band 
structure of (Ga,Mn)As only in bulk systems . 5161718 Some calculations included the ef- 
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fects of biaxial strain, spin-orbit coupling, and exchange correlation in a parabolic band 
approximation . 9 ! 10 ! 11 In the case of quantum wells, (Ga,Mn)As multilayers and superlat- 
tices, a self-consistent calculation has been performed for parabolic h eavy holes subbands.— 
Self-consistent calculations have also been performed in Refs. liol and lllf assuming isotropic 
effective masses. To the exception of Ref. [12|, where a Monte Carlo simulation is used, the 
electronic structure calculation assumes a homogeneous density of magnetic moment, as well 
as homogeneous negative charge concentration (due to the ionized Mn atoms). Relaxing this 
approximation implies in considering a multiple-scattering treatment, what is outside the 
scope of the present work. The homogeneous approximation for the density of magnetic and 
Coulomb scattering centers (localized magnetic moments and ionized impurities) provides 
important information concerning the carriers charge and spin distributions. 

Here we present a self-consistent LK k-p calculation for GaMnAs/GaAs multilayers and 
superlattices. As described below, we adopt a super-cell calculation which is an extension of 
the LK method to treat the cases of quantum wells and superlattices (SL). The structure we 
consider consists of substitutional Mn ions uniformly distributed in Gao.95Mno.05As layers 
of width d\, with a hole concentration equivalent, in bulk, to lxlO 20 cm -3 , assumed to be 
metallic and ferromagnetic, at T=0K. The DMS are separated by non-magnetic GaAs layers 
of width c?2- Before the first DMS layer and after the last one, GaAs lateral layers of width 
s complete the structure, as shown in Fig. ^ The super-cell model consists of placing this 
"active part" , described so far, between thick layers of a large gap material, assumed here 
to be AlAs, and treating the whole system as a superlattice. 

The final spin and charge configuration contains a complete information about the com- 
position of heavy holes and light holes in the structure. Carriers are anti-parallel heavy holes, 
and anti-parallel light holes, in a lesser amount. A strong spin polarization is observed, and 
the charge is concentrated mostly on the GaMnAs layers. The charge and spin distributions 
are analyzed in terms of their dependence on the number of multilayers, the widths d\ and 
<i 2 , as well as the widths of the lateral GaAs layers. 
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II. MAGNETIC INTERACTION IN THE CONFINED LK MODEL 



The interaction between free holes and the localized magnetic moments is well described 
by the Kondo-like term 

V mag (r) = -/]Ts(r) • S(R,)5(f- R 4 ), (1) 

i 

where / is the sp — d interaction. The localized spin of the Mn ion Si at position Ri is 
treated as a classical variable, since it results of the five 3d electrons obeying Hund's rule, 
and no hybridization with carriers is considered, due to the high difference in energies. s(r) 
is the spin operator of the carrier at position r. At zero temperature, assuming a complete 
alignment of the localized magnetic moment, i.e., S(R i ) = S, we have: 

V mag {f) = -IS ■ s(f) £ 5(f - Ri) = —S ■ a Pl {r). (2) 

i 

In the last expression we used a = ia x + jo y + ka x z to denote the three Pauli matrices; 
Pi(r) is the density of magnetic impurities. Assuming a homogeneous distribution of the 
localized magnetic dipoles inside the DMS layers , we have Pi(r) ~ xN g(z), where N is 
the the density of cations, x is the substitutional concentration of Mn, and g(z) — 1 if z lies 
inside a DMS layer, g(z) = otherwise. In that case the magnetic interaction becomes: 

V mag (z) = ~N (3g(z)M-a. (3) 

In Eq.fjSJ) we have explicitly taken into consideration that carriers are holes, replacing / by 
p. For GaMnAs, N Q (3 = -1.2 eV, M For electrons, N a = 0.2 eV.^ Were the spin of the 
particles well defined, this term would represent, in bulk, a shift on the top (bottom) of 
the valence (conduction) band. This is not the discussed above, for light holes and 

split-off holes in GaAs. 

It is well known that the valence bands in GaAs at the T point split, due to spin-orbit 
interaction, into four j = 3/2 states belonging to the Tg representation, and two j = 1/2 
states belonging to the representation. These two are separated from the r 8 states by 
the spin-orbit energy A, which is 340 meV in GaAs.-^ We adopt the notation | j, rrij > to 
represent the r 8 and the T 7 states, making use of the fact that these states are also eigenstates 
of the total angular momentum operator J = L + s , with eigenvalue j, and simultaneously 
eigenstates of its z-component, J z , corresponding to the eigenvalue rrij. Notice that the 
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Tg states | vl > and | v2 > both having j = 3/2 but with rrij = 3/2 and rrij = —3/2, 
respectively, are called heavy holes states, having their spins well defined, being "up", i.e., 
aligned with J, and "down", anti-aligned. They can be represented by using as a basis the 
three p-type states in the directions x, y, and z:— 

\vl> = \\\>=^\ x + l V» ( 4 ) 

\v2> = \ -- > =±=(\ X -iy[ > ) (5) 

The other states, the light hole states | v3 > and | t>4 > corresponding also to j = 3/2, but 
with rrij = 1/2 and rrij = —1/2, and the split-off states | v5 > and | v6 > corresponding to 
j = 1/2 with rrij = 1/2 and rrij 

| v3 > = 

| vA > = 

| v 5 > = 

| v6 > — 

The kinetic, Hartree, and the exchange-correlation terms appear in this formalism as the 
well known components of the LK matrix. 4 In multilayers and superlattices the mismatches 
of the valence and conduction bands, which play the roles of confining potentials, are to be 
added. Differences on the lattice parameters introduce the additional terms of the strain. 
Here we have also to introduce the magnetic potential given by Eq.fjHI). Two distinct cases 
may be considered, since a break in the T^ symmetry occurs: magnetization "in-plane" 
(occurring parallel to the interfaces) Mil = M x i + M y j, and "perpendicular to the plane", 
= M z k. The latter is assumed to be in the growth direction, here considered as the z- 
axis. Notice that in the presence of a confining potential created by the interfaces, V c (z), the 
operators J x and J y no more commute with the LK Hamiltonian. In bulk, however, within 
the homogeneous magnetization approach, there is no distinction between these two cases. 
Once we are not trying to explain the origin of the ferromagnetic order in these systems, the 
occurrence of the magnetization being "in plane" or "perpendicular-to-the-plane" is assumed 
to be provided by an external weak magnetic field, which does not interfere on the electronic 
structure, directly. 
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= —1/2, do not have well defined spins: 

I 22 >= 7^ -2 \ *> + \x + W l>) (6) 

| ^I>=-L(|x-n/T>+2Ui>) (7) 

\\\>=^{\z]>^\x + iy[>) (8) 

l M >= 73 (-|a;- ^ T> + |2;i>) (9) 



For the sake of obtaining the LK matrix for heterostructures, it is necessary to calculate 
first the matrix elements < j, m\V mag \j', m' > for each constituent DMS layer, in bulk. This 
is easily done by observing that 



M ± -a\ t> = M z \ |>, 
M ± -a\ i> = -M z \ |>, 



(10) 
(11) 



and 



My ■ a\ t> = (Mj + iMy)\ [>= M + | |>, 
M v a\ [> = (M x -iM y )\ t>=M_| T>, 



(12) 
(13) 



Making use of the approximation given in Eq.® and the results in Eqs. (jl(JH13|) . we have 
for the LK matrix of V mag at the T point: 



V, 



mag 
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(14) 



During the last years the LK model has been adapted to quantum wells and superlattices 



(SL), as described in Refs. 
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18U19L We adopt that approach, using a super-cell model. This 
means that we consider an unit cell consisting of the active region plus a thick insulator layer. 
The number of DMS layers in the unit cell can be varied at will. We assume, then, an infinite 
SL in the [001] direction. The multiband effective-mass equation (EME) is represented with 
respect to plane waves with wavevectors K = {2n/a)l (I an integer and a the SL period) e qual 
to the reciprocal SL vectors. A detailed description of the method can be found in Ref. |22t 
The rows and columns of the 6x6 LK Hamiltonian relate to the Bloch-type eigenfunctions 
\j,rrij, k > of the T 8 heavy- hole bands, and the spin-orbit-hole band; k denotes a vector 
of the first Brillouin zone. Expanding the EME with respect to plane waves < z \ K > 
means representing this equation in terms of the Bloch functions < x | j,rrij,k + Ke z >. 
For a Bloch function < z \ E, k > of the SL corresponding to energy E and wavevector k, 
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the EME takes the form: 

< J, rrij, k,K\T + H s + V het + V C + V xc + V mag \ f, mj, k,K' > x 

<j',m' p k,K' \E,k>= E{k) < j, rrij, k, K \E,k>, (15) 

where T is the unperturbed kinetic energy term generalized for a heterostructure,— H$ is 
the strain energy term originating from the lattice mismatch, Vhet is the square potential due 
to the difference between energy gaps, V xc is the exchange-correlation potential, and Vc is 
the sum of the Hartree potential with the ionized acceptor potential. Finally, V mag is given 
by Eq.fjHI), for each material. The Luttinger parameters and the other terms appearing in 
the secular equation are to be taken for each epitaxial layer of the SL.— For instance, in the 
case of the magnetic interaction we have: 

< j, m, k, K\V mag \f, m', k', K' >=< k, K\V^f m 'g{K' - K)\k> , K> >, (16) 

where the integral 

g(K'-K) = - f" e- iKz g{z)e iK ' z (17) 
a Jo 

is performed in a DMS layer of width d. 

The self-consistent potentials and the charge densities are obtained by solving the multi- 
band EME equation and the Poisson equation: 

47re 2 1 

< j, rrij, k,K\V c \ f, m'j, k, K' >= — — - < K \ p + + p \ K' > 5^,8^^ (18) 

J K \K — K \ 3 

where k is the dielectric constant of the host, and p + and p~ are the density of charge of 

holes and acceptors, respectively, expressed in plane-wave representation. 



III. RESULTS 

The DMS layers work effectively as barriers or wells for spins parallel (up) and antipar- 
allel (down) to the local average magnetization, depending on the sign of iVo/3 for valence 
band, and N a for conduction band. These DMS layers are assumed to be ferromagnetic 
and metallic, with a 3-D equivalent hole density p — 1 x 10 20 cm -3 , a substitutional Mn 
concentration of 5%, at temperature T= K, and an average magnetization (M)=5/2. 

In Fig. |2]we present (a) the valence band structure (hole binding energy) and (b) potential 
profiles for a system consisting of six DMS layers, with 0^=20 A and d,2=30 A. Energies are 
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reckoned from the top of the Coulomb barrier, as in Ref. |22|. The magnetization is assumed 
to be in the z-direction. The subbands are hybrid states, since they are a mixing of all kinds 
of holes. However, at the T-point the lowest lying states have a dominant component. For 
instance, the first three states are almost entirely heavy holes "down". The mixing becomes 
stronger in the more excited states, and as we go out of the T-point. In Fig. |2]we named 
the band by its dominant component at the T-point. Here, "up" and "down" refer to the 
sign of rrij, the z-component of the total angular momentum. In other words, "up" means 
parallel to the average magnetization, while "down" means anti-parallel. The T-A (r~Z)- 
line corresponds to wave vectors k perpendicular (parallel) to the SL axis. The Fermi energy 
is also indicated. Strong non-parabolicity arises in the subbands along the (r-A)-line, which 
leads to remarkable anti-crossing behavior. We also observe that several levels are occupied, 
most of them are heavy and light hole down. It is possible to understand this behavior 
by observing Fig. EJb), which shows the self-consistent hole band potential profile for each 
carrier. The confinement for heavy and light holes down are deeper than the cases where 
the z-components of the total angular momentum are up. 

Fig. El shows the carriers distribution for the structure with six DMS layers described 
in Fig. |21 The results are shown for lateral GaAs layers of width s =0, 10, 20, 30, 40, 
50 and 60 A. We notice that, except for the heavy hole with rrij = 3/2 (up), carriers 
concentrate on the DMS layers. The heavy holes up concentrate in the non-magnetic region 
as a consequence of the strong magnetic repulsion. After s = 50A the carriers distribution 
becomes independent of the lateral GaAs widths. 

Fig. 0] shows the carriers density for active regions consisting of 1 to 9 DMS layers. We 
fixed <ii=20 A and g?2=30 A. The charge density is plotted for each component. As before, 
"up" and "down" refer to the sign of the rrij component of the total angular momentum. 
These results are consistent with the electronic structure shown in Fig. El since the lowest 
levels are mostly heavy and light holes down. Observe that in all multilayered structures (2 
to 9 DMS layers) the charge is concentrated almost entirely inside the DMS layers. Again, 
the exception is the heavy hole up component, whose density is higher in the non-magnetic 
regions for the reasons explained above. 

Up to now the distribution of charge in the multilayered structure has been shown in 
terms of carriers total angular momentum components. For the sake of completeness, it is 
also interesting to know how the charge is distributed in terms of the spin polarization. This 
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can be easily obtained, since the light holes are mixed states of up and down spins with 
defined probabilities. The results for these multilayers are shown in FigJS] in terms of the 
spin components of the charge density There is a strong spin polarization, dominated by 
heavy and light holes down spins. 

In Fig. El we analyze the dependence of the effective two-dimensional hole concentration, 
N2D, which is the integrated 3-D density on the z-direction for each component. The DMS 
layers widths are changed, while keeping fixed the width of the GaAs layers, d,2=30 A. The 
calculation is performed for 6 DMS layers. For the heavy and light hole down, the 2-D 
densities are much higher then those associated to positive values of rrij. The difference 
is one order of magnitude, in the case of heavy holes, due to its higher effective mass and 
profile potential that provide a higher occupation inside the DMS layers. 

Finally, in Fig. [7|we show N2D as a function of the non-magnetic layer width, d,2, with 
di=20 A, for each carrier, heavy and light hole (down and up). The calculation is performed 
for 4, 5, 6 and 7 DMS layers. 



IV. CONCLUSIONS 



In summary, we have investigated the electronic structure of Gao.95Mno.05As/GaAs mul- 
tilayers by using a super-cell model in the framework of the Luttinger-Kohn k.p approxi- 
mation. The unit in the super-cell is a AlAs/n(Gao.95Mno.o5As/GaAs)/AlAs structure, with 
n being the number of DMS layers, grown in the [100] direction. The DMS layers are as- 
sumed to be metallic and ferromagnetic, at T=0K. A small magnetic field is applied in the 
growth direction to guarantee the magnetization to be perpendicular to the plane, but it 
does not affect the electronic structure. Several subbands are occupied. They are mostly 
heavy holes with rrij = —3/2, the carrier density being higher in the DMS layers. However, 
a non- negligible part of carriers are light holes with rrij = —1/2, also concentrated on the 
same region as the heavy holes. Actually, what is observed is the appearance of two differ- 
ent channels, one with the z-component of the total angular momentum aligned with the 
average magnetization, which is highly concentrated on the DMS layers, and another, with 
opposite orientation whose distribution is more concentrated on the non-magnetic layers of 
the active region, the latter channel being much less dense. These results are in qualita- 
tive agreement with those obtained in Ref. showing the electronic structure calculations 
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in digital ferromagnetic heterostructures of GaAs/(Ga,Mn)As and AlAs/(Ga,Mn)As/GaAs 
using density-functional theory in the local spin-density approximation. As a consequence 
of both magnetic and Coulomb interactions - here included the hole-hole interaction - the 
spin-polarized charge tends to be slightly non-periodically distributed as the number of DMS 
layers increases. Carriers polarized anti-parallel to the magnetization tend to concentrate a 
little more at the borders of the structure, while carriers polarized in the opposite direction 
tend to concentrate on the middle. 

From the point of view of the observation of a ferromagnetic multilayer with a high 
transition temperature (we recall that we have assumed in our calculation T=0K), it is 
interesting to have a strong "spin-polarized" charge density in the DMS layers, while keeping 
some charge in between to guarantee the inter-layer interaction . 24 ' 25 However, for obtaining 
a high mobility spin-polarized current, a structure should be grown in which the in-plane 
transport is realized with spin-polarized carriers concentrated in a region of high mobility, 
away from the scatterers. Therefore, the ideal distribution of charge and spin depends on 
the purpose of the structure. The results we have obtained points to the possibility of 
engineering the spin-polarized charge distribution by the right choice of the magnetic layers 
and the band mismatches with the non-magnetic spacers. 
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FIG. 1: The model structure: the DMS layers of width dl are separated by non-magnetic GaAs 
layers of width d 2 . At the left of the first DMS layer and at the right of the last one, GaAs lateral 
layers of width s complete the "active" parts. The "active" parts are separated, themselves, by 
thick AlAs barriers. 

FIG. 2: (a) Valence band structures for GaAs/GaAs:Mn multiple QWs along high symmetry lines 
T-Z and T-A, with d\ = 20 A and d 2 = 30 A. We named the band by its dominant component 
at the T-point. Solid lines correspond to subbands which, at the T-point, are mostly heavy holes 
down, dashed lines to light hole down, dotted lines to heavy hole up, and dotted-dashed lines 
to light holes up. The Fermi energy is also indicated (short-dashed line), (b) the corresponding 
potential profile for each carrier, down and up. Energies are reckoned from the top of the Coulomb 
barrier. 

FIG. 4: Particle density of heavy and light holes for a structure consisting from 1 to 9 DMS layers. 
Widths d\ and d 2 as in Fig^ 

FIG. 5: Same as in Fig. |31 for each carrier spin component. 



FIG. 6: N2D as a function of the DMS layer width d\. The width of the non-magnetic layers are 
fixed, d 2 = 30 A. 



FIG. 7: N2D as a function of the non-magnetic layer width d 2 , for 4, 5, 6 and 7 DMS layers. The 
widths of the DMS layers are fixed, d x = 20 A. 

FIG. 3: Density distribution of heavy and light holes for six DMS layers with di = 20A, d 2 = 30A, 
and lateral GaAs widths s = 0, 10, 20, 30, 40, 50 and 60A. 
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